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We compare results of the fc t -factorization approach and the next-to-leading order collinear- 
factorization approach for photon-jet correlations in pp and pp collisions at RHIC and Tevatron 
energies. We discuss correlations in the azimuthal angle as well as in the two-dimensional space 
of transverse momentum of photon and jet. Different unintegrated parton distributions (UPDF) 
are included in the fct-factorization approach. The results depend on UPDFs used. The standard 
collinear approach gives cross section comparable to the fct-factorization approach. For correlations 
of the photon and any jet the NLO contributions dominate at relatively small azimuthal angles as 
well as for asymmetric transverse momenta. For correlations of the photon with the leading jet 
(the one having the biggest transverse momentum) the NLO approach gives zero contribution at 
<j>- < tv/2 which opens a possibility to study higher-order terms and/or UPDFs in this region. 

PACS numbers: 12.38.Bx, 13.60.Hb, 13.85. Qk 



I. INTRODUCTION 

The jet-jet correlations are interesting probe of QCD 
dynamics p}. Recent studies of hadron-hadron correla- 
tions at RHIC 3 open a new possibility to study the 
dynamics of jet and particle production. The hadron- 
hadron correlations involve both jet-jet correlations as 
well as complicated jet structure. Recently a preliminary 
data on photon-hadron azimuthal correlation in nuclear 
collisions were also presented |3j. In principle, such cor- 
relations should be easier for theoretical description as 
here only one jet enters, at least in leading order pQCD. 
On the experimental side, such measurements are more 
difficult due to much reduced statistics as compared to 
the dijet studies. 

Up to now no theoretical calculation for photon-jet 
were presented in the literature, even for elementary col- 
lisions. In leading-order collinear-factorization approach 
the photon and the associated jet are produced back- 
to-back. If transverse momenta of partons entering the 
hard process are included, the transverse momenta of 
the photon and the jet are no longer balanced and finite 
(non-zero) correlations in a broad range of relative az- 
imuthal angle and/or in lengths of transverse momenta 
of the photon and the jet are obtained. The finite cor- 



relations can be also obtained in higher-order collinear- 
factorization approach (4)]. According to our knowledge 
no detailed studies for present accelerators were pre- 
sented in the literature. 

In contrast to the coincidence studies the inclusive dis- 
tributions of photons were studied carefully in pQCD 
up to the next-to-leading order @. Similar studies were 
crformed recently also in the fct-factorization approach 
0] ■ A rather good description of direct-photon inclu- 
sive cross sections can be obtained in both approaches. 
The fc t -factorization approach offers a relatively easy 
method to calculate the photon-jet correlations p. 

The fct-factorization approach was used recently to sev- 
eral high-energy reactions, including heavy quark pair 
photo- [1, Hand hadroproductio n flO . , charmonium 
production [l2l [l3j , inclusive Zq [14| and Higgs [lH, flU ] 
production. 

In the present paper we shall compare results ob- 
tained in the leading-order k t -factorization approach 
and the ncxt-to-lcading order collinear-factorization ap- 
proach. We shall discuss which approach is more ade- 
quate for the different regions of phase space. We shall 
present corresponding results for proton-proton scatter- 
ing at RHIC and proton-antiproton scattering at Teva- 
tron. 
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II. FORMALISM 
A. 2^2 contributions with unintegrated parton distributions 

It is known that at midrapidities and at relatively small transverse momenta the photon-jet production is dominated 
by (sub)processes initiated by gluons. In Fig[T]we show basic diagrams which appear in the fct-factorization approach 
to photon-jet correlations. 
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FIG. 1: Basic diagrams of fc t -factorization approach to photon-jet correlations. 



In the /cj-factorization approach the cross section for a simultaneous production of a photon and an associated jet 
in the collisions of two hadrons (pp or pp) can be written as 

da hlh2 ^ lk f cfki.t d 2 kx.t 1 



, dyxdyi ; ~Ta~2? ^l^i-^kl 2 

■ S 2 (k ltt + k 2:t -pi,t -p2,t)^t(xi,kl t ,fi 2 1 )J 7 3 (x 2 ,k^ t ,nl) , (1) 

where Fi(xi, k\ t , fj, 2 ) and Tj(x 2 -, k\ t , //|) are so-called unintegrated parton distributions. The longitudinal momentum 
fractions are evaluated as 

xx = (mi t e +Vl + m 2 te +y2 ) / \fs , 

x 2 = (m u e- yi + m 2t e- y2 )/^. (2) 

We shall return to the choice of the factorization scale in the next section. Its role is completely different in different 
approaches i.e. different choices of UPDFs. A special attention will be devoted to the Kwiecihski UPDF and the role 
of the scale paramater. 

If one makes the following replacement 

Ti(xi,k\ t ) -> xipi{xi)5(kl t ) , 

Tj{x 2 ,k\t) -> x2Pj{x2)5{kl t ) (3) 

then one recovers the standard leading-order collinear formula. 

The final partonic state is jk — 75, jq. The matrix elements for corresponding processes are discussed in Appendix 

A. 

The inclusive invariant cross section for direct photon production can be written as 

da hl h 2 ^ _ f, d 2 k ht d 2 k 2 , t ( , . 

dyid 2 p^ t ~ J V2 7T 7T [ - )l P2, t =K t +k2, t -Pi, t W 

and analogously the cross section for the associated parton (jet) can be written as 

d<7 hlh ^ k r d 2 k lit d 2 k 2 .t , )} ft ., 

dyid 2 Pht ~ J Vl 7T ir [ - )l pi, t =ki, t +k2, t - P \ t ■ W 
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Let us return to the coincidence cross section. The integration with the Dirac delta function in Eq.JT]) 



dy\dy 2 



-(•■•)<5 2 (-) 



7T 7T 

can be performed by introducing the following new auxiliary variables: 

Qt = ki.t + k 2 . t , 
Qt = - k 2 . t ■ 

Then our initial cross section can be written as: 

d(7h\h-2 — >7 ,parton 



If 1 

4 / dy x dy 2 ^dq^d<j) qt (...)\^ t= p t 



d 2 Pi tt d 2 p 2 ,t 4 

Above P t = pij +P2.t- The factor i before integrand on the rhs comes from the Jacobian of the (&i,t, k 2 



transformation (see [6|) 



(6) 



(7) 

(8) 
{Qt,qt) 



B. 2 — > 3 contributions in NLO collinear-factorization approach 

Up to now we have concentrated only on processes with two explicit hard partons {"fk) in the fcf-factorization 
approach. It is of interest to compare the results of our approach with those of the standard collinear next-to-leading 
order approach. In this section we discuss processes with three explicit hard partons. In Fig|5]we show diagrams 
for 2 — > 3 subprocesses included in our calculations. In the following we assume particle No.l to be a photon. Then 
particle No. 2 is g,q and q, depending on the subprocess. 
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FIG. 2: Diagrams for NLO collinear-factorization approach for photon-jet-jet production. 

The cross section for h\h 2 — > jklX processes can be calculated according to the standard parton model formula 

dcr h ih 2 ^jki = X! [dx 1 dx 2 Pi{xi,fj/ i )pj(x2,lJ?)da i j^ryki. (9) 
ijk 

The elementary cross section can be written as 

d&ij-^kl = T^Mij^ki^dRz , (10) 
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where the three-body phase space element reads 



dR 3 = {^) i 5 i (Pa+p b -jzp i )f{ jy:^ ■ (11) 

i=l i=l /£/ ^ /7r J 



This element can be expressed in an equivalent way in terms of parton rapidities 

dR 3 = (2,f8\p a + Pb - J» f[ ■ (12) 

f~[ fJi (47r)(27r) 2 

The last formula is useful for practical applications. Now the cross section for hadronic collisions can be written in 
terms of 2 — > 3 matrix element as 

do = ^2dy 1 d 2 p 1 ,tdy 2 d 2 p2 : tdy3 , , 3 \ ^x 1 p i (x 1 , ^ 2 )x 2 Pj{x 2 , u 2 )\M io -, lk i\ 2 , (13) 

ijkl \ \ s 

where the longitudinal momentum fractions are evaluated as 

3 



4=E^e"^. (W) 



3 
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Repeating similar steps as for 2-^2 processes we get finally 

do- = ^-^^x lPl {x 1 ,^ 2 )x 2 p ] (x2,^ 2 )\M lj ^ki\ 2 Pi,tdpi,tP2,tdp2,td(f>-dy 1 dy2dy 3 , (15) 

ijkl 

where the relative azimuthal angle between the photon and the associated jet (4>~) is restricted to the interval (0,7r). 
The last formula is very useful in calculating the cross section for particle 1 and particle 2 correlations. 



III. RESULTS 



In this section we shall present results for RHIC and Tevatron energies. We use UPDFs from the literature. There 
are only two complete sets of UPDF in the literature which include not only the gluon distributions but also the 
distributions of quarks and antiquarks: 

(a) Kwiecihski [17}, 

(b) Kimber-Martin-Ryskin [iH |. 

For comparison we shall include also the unintegrated parton distributions obtained from the collinear ones by the 
Gaussian smearing procedure. Such a procedure is often used in the context of direct photons [19l . [20(. Compar- 
ing results obtained with those Gaussian distributions and the results obtained with the Kwiecihski distributions 
with nonperturbative Gaussian form factors will allow to quantify the effect of UPDF evolution as contained in the 
Kwiecihski evolution equations. What is the hard scale for our process? In our case the best candidate for the scale 
is the photon and/or jet transverse momentum. Since we are interested in rather small transverse momenta the 
evolution length is not too large and the deviations from initial fc t -distributions (assumed here to be Gaussian) should 
not be too big. 

At high energies one enters into a small-x region, i.e. the region of a specific dynamics of the QCD emissions. In 
this region only unintegrated distributions of gluons exist in the literature. In our case the dominant contributions 
come from QCD-Compton gluon — quark or quark — gluon initiated hard subprocesses. This means that we need 
unintegrated distributions of both gluons and quarks/antiquarks. In this case we take such UGDFs from the literature 
and supplement them by the Gaussian distributions of quarks/antiquarks. 

Let us start from presenting our results on the (pi,t,P2,t) plane. In Fig[3]we show the maps for different UPDFs 
used in the fct-factorization approach as well as for NLO collinear-factorization approach for pi t t,P2,t € (5, 20) GeV 
and at the Tevatron energy W = 1960 GeV. In the case of the Kwiecihski distribution we have taken bo = 1 GeV -1 
for the exponential nonperturbative form factor and the scale parameter /i 2 = 100 GeV 2 . Rather similar distributions 
are obtained for different UPDFs. The distribution obtained in the NLO approach differs qualitatively from those 
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obtained in the £; t -factorization approach. First of all, one can see a sharp ridge along the diagonal px,t = P2.t- This 
ridge corresponds to a soft singularity when the unobserved parton has very small transverse momentum p 3 t . As will 
be clear in a moment this corresponds to the azimuthal angle between the photon and the jet being (/>_ = tt. Obviously 
this is a region which cannot be reliably calculated in collinear pQCD. There are different practical possibilities to 
exclude this region from the calculations. The most primitive way (possible only in theoretical calculations) is to 
impose a lower cut on transverse momentum of the unobserved parton p^ t t- Secondly the standard collinear NLO 
approach generates much bigger cross section at configurations asymmetric in p\ tt and P2,t- We shall return to this 
observation in the course of this paper. 





FIG. 3: Transverse momentum distributions da/dpi t tdp2,t at W = 1960 GeV and for different UPDFs in the fct-factorization 
approach for Kwiecihski (bo — lGeV -1 , /i 2 = lOOGeV 2 ) (a), BFKL (b), KL (c) and NLO 2-^3 collinear-factorization approach 
including diagrams from Figf2](d). The integration over rapidities from the interval -5 < J/i, yi < 5 is performed. 



As discussed in Ref.[6| the Kwiecihski distributions are very useful to treat both the nonperturbative (intrinsic 
nonperturbative transverse momenta) and the perturbative (QCD broadening due to parton emission) effects on the 
same footing. In Fig0] we show the effect of the scale evolution of the Kwiecihski UPDFs on the azimuthal angle 
correlations between the photon and the associated jet. We show results for different initial conditions (bo = 0.5, 1.0, 
2.0 GeV -1 ). At the initial scale (fixed here as in the original GRV 2l| to be /i 2 = 0.25 GeV 2 ) there is a sizeable 
difference of the results for different bo- The difference becomes less and less pronounced when the scale increases. At 
= 100 GeV 2 the differences practically disappear. This is due to the fact that the QCD-evolution broadening of the 
initial parton transverse momentum distribution is much bigger than the typical initial nonperturbative transverse 
momentum scale. 

In FigjS] we show corresponding azimuthal angular correlations. In this case integration is made over transverse 
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FIG. 4: (Color on line) Azimuthal angle correlation functions at (a) RHIC, (b) Tevatron energies for different scales and 
different values of bo of the Kwieciriski distributions. The solid line is for bo = 0.5 GeV -1 , the dashed line is for bo = 1 GeV - 
and the dotted line is for bo = 2 GeV -1 . Three different values of the scale parameters are shown: /i 
(the bigger the scale the bigger the decorellation effect, different colors on line). In this calculation pi tl 
J/i,2/2 G (-5,5). 



= 0.25, 10, 100 GeV 2 
P2,t e (5,20) GeV and 



momenta p\,uP2,t & (5,20) GeV and rapidities 2/1,2/2 G (—5,5). 
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FIG. 5: Photon-jet angular azimuthal correlations da/d(j>- for proton-antiproton collision at W = 1960 GeV for different 
UPDFs in the fct-factorization approach for the Kwieciriski (solid), BFKL (dashed), KL (dotted) UPDFs/UGDFs and for the 
NLO collinear-factorization approach (thick dashed). Here 2/1,2/2 £ (—5,5). 



The singularity in NLO pQCD at </>_ = it is strongly correlated with the sharp ridge in FigJ3]d. This is demonstrated 
in Fig[S] where we present the results of azimuthal correlation function obtained for different cuts on ps y t- The cut 
modifies only the region of relative azimuthal angles close to it. We wish to stress in this context that there are no 
singularities of the ridge type in the k t -factorization approach. 

The small transverse momenta of the unobserved jet contribute to the sharp ridge along the diagonal pij = P2.t- 
It is therefore difficult to distinguish these three-parton states from the states with two partons. The ridge can 
be eliminated in calculation by imposing a cut on the transverse momentum of the third (unobserved) parton. In 
experiments there is no possibility to impose such cuts and other methods must be used. We shall return to this point 
later in this paper. 

In FigJT] we show angular azimuthal correlations for different relations between transverse momenta of outgoing 
photon and partons: (a) with no constraints on p%j, (b) the case where p2,t > P3,t condition (called leading jet 
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FIG. 6: Photon-jet angular azimuthal correlations dcr/dcj)- for proton-antiproton collision at W = 1960 GeV for the NLO 
collinear-factorization approach and different cuts on pz,t- Here 2/1,2/2 6 (—5,5). 



condition in the following) is imposed, (c) p2,t > P3,t and an additional condition pi^ > P3,t- The results depend 
significantly on the scenario as can be seen from the figure. 
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FIG. 7: Angular azimuthal correlations for different cuts on the transverse momentum of third (unobserved) parton in the 
NLO collinear-factorization approach without any extra constraints (dashed), p3, t < p2,t (solid), p3, t < p2,t and j>3, t < pi tt in 
addition (dotted). Here 3/1,3/2 € (—5,5). 



In FigJ5]we show transverse momentum distribution da / dp\ t tdp2,t for the same extra conditions imposed before in 
Figj7]for azimuthal angle correlations. Imposing the condition that the associated jet is the leading jet (j>2,t > P3,t) 
causes that the large part of the phase space p\ t t < ^P2,t is not available in the next-to-leading approach. If one 
imposes in addition that pi^ (photon) > p^j (unobserved jet) then also the p2,t < %Pi,t region becomes excluded for 
the NLO approach. These NLO-excluded regions are therefore regions sensitive to higher-order corrections in pQCD. 

In general, the correlations between the photon and the jet depend strongly on all kinematical variables - transverse 
momenta, azimuthal angles, etc. In order to expose this better, in FigJ5]we define windows in the (j>i,t,P2,t) plane 
which will be used in the following to study the azimuthal correlations. At lower energies (as for RHIC) a region 
of rather low transverse momenta is more adequate (left figure). At larger energies (as for Tevatron) also region of 
somewhat larger transverse momenta can be of interest (right figure). The notation shown in the figure will be used 
for brevity in the rest of this paper for easy reference. 

In FigfTniwe show angular azimuthal correlations dcr/dcj)- at RHIC energy — 200 GeV for the Kwiecihski UPDFs 
in the fc t -factorization approach with on-shell and off-shell matrix elements and for the NLO collinear-factorization 
approach with extra leading jet condition p^ tt < P2,t- Here transverse momentum of the photon and that of the 
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(c) 

FIG. 8: Transverse momenta distribution do j dp\, t dpi,t for different constraints on the transverse momentum of third (un- 
observed) parton in the NLO collinear-factorization approach with no constraints on pz,t (a), pz,t < P2.t (b), pz t t < P2,t and 
P3,t < Pi,t (c). All 2 — » 3 processes shown in Fig[2]were included. Here 2/1,2/2 £ (—5,5). 



associated jet pi } t,P2,t belong to the interval (5, 20) GeV. There is almost no difference between results obtained with 
off-shell and on-shell (see Appendix A) matrix elements. 

In Fig[Tl]we show analogous angular distributions as in FigfTUlbut for Tevatron energy y/s = 1960 GeV. In Fig[T2l 
we show angular correlations for a restricted range of rapidities \yi\, \y%\ < 0.9 (corresponding to the present Tevatron 
apparatus) of the photon and the correlated jet. Limiting to midrapidities does not change the shape of azimuthal 
correlations significantly. 

In Fig |13l we show similar distributions as in Fig llll but for transverse-momentum windows spanned over broader 
range of transverse momenta Pi,t,P2,t € (20, 80) GeV for the photon and the jet. We observe slightly faster decrease 
of the £; t -factorization cross sections for larger p\ it and pi.t- 

The standard collinear approach can be applied only in the region which is free of singularities. In order to eliminate 
the regions where the pQCD calculation is not reliable some cuts on the measured transverse momenta must be applied. 
The simplest method is to use cuts shown in Fig |141 Mathematically this means that pu > Pcuti P2,t > Pent and 

\p l ,t-p 2 ,t\ > A s . (16) 

We shall call the last cut a scalar cut for further easy reference. In Fig[15] we show azimuthal angle correlation 
function for different values of the scalar cut A5 = 0, 1,2,3 GeV. Clearly the NLO singularity at <fi- = ir can be 
removed by imposing the cut. However, the cut lowers also the /c t -factorization cross section. 
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FIG. 9: The definition of the windows in (pi,t,P2,t) plane for RHIC energy y/s = 200 GeV (a) and for Tevatron energy 
yfs = 1960 GeV (b). 
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FIG. 10: Angular azimuthal correlations dcr/dcj)- at y/s = 200 GeV for Kwiecihski on-shell ME (solid), Kwiecihski off-shell 
ME (thick dotted), NLO collinear with no cuts on p3,t (dashed) and NLO collinear with cut on p3 it < p2,t (dash-dotted). Here 
Vi,V2 € (-5,5). 



We have also tried another option to cut off the singularity: 

\Pl,t +P2,t \ > Ay 



(17) 



This type of the cut will be called vector one for brevity. In Fig llGl we show corresponding photon-jet azimuthal angle 
correlation function with different values of the cut Ay = 0, 1, 2, 3 GeV. The situation here is very similar to that for 
the scalar cut. 



IV. CONCLUSIONS 



tion observables in proton-proton (RHIC) and proton- 



We have performed for the first time the lacking 
in the literature calculation of the photon-jet correla- 
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FIG. 11: Angular azimuthal correlations dcr/dcj)- at yfs — 1960 GeV for Kwieciriski on-shell ME (solid), Kwieciriski off-shell 
ME (thick dotted), NLO collinear with no cuts on p3,t (dashed) and NLO collinear with cut on p3. t < p2,t (dash-dotted). Here 
Vl,V2 G (-5,5). 
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FIG. 12: Angular azimuthal correlations dcr/dcj)- at y/s = 1960 GeV for Kwieciriski on-shell ME (solid), Kwiecihski off-shell 
ME (thick dotted), NLO collinear with no cuts on p^,t (dashed) and NLO collinear with cut on ps >t < p2,t (dash-dotted). Here 
1 2/1 1 , | J/2 j < 0.9. 



antiproton (Tevatron) collisions. Up to now such corre- have concentrated on the region of small transverse mo- 
lations have not been studied experimentally either. We menta (semi-hard region) where the k t -factorization ap- 
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FIG. 13: Angular azimuthal correlations da/dfi- at y/s — 1960 GeV for Kwieciriski on-shell ME (solid), Kwieciriski off-shell 
ME (thick dotted) NLO collinear with no cuts on p3,t (dashed) and, NLO collinear with cut on p3, t < p2,t (dash-dotted). Here 
2/1,1/2 € (-5,5). 
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FIG. 14: Diagram showing excluded region in (pi,t,P2,t) plane. 



proach seems to be the most efficient and theoretically 
justified tool. We have calculated correlation observables 
for different unintegrated parton distributions from the 
literature. Our previous analysis of inclusive spectra of 
direct photons suggests that the Kwiecihski distributions 
give the best description at low and intermediate ener- 
gies. We have discussed the role of the evolution scale of 
the Kwiecihski UPDFs on the azimuthal correlations. In 
general, the bigger the scale the bigger decorrelation in 



azimuth is observed. When the scale fj 2 ~ pi (photon) 
~ pi (associated jet) (for the kinematics chosen /i 2 ~ 
100 GeV 2 ) is assumed, much bigger decorrelations can 
be observed than from the standard Gaussian smearing 
prescription often used in phenomenological studies. 

The correlation function depends strongly on whether 
it is the correlation of the photon and any jet or the cor- 
relation of the photon and the leading-jet which is consid- 
ered. In the last case there are regions in azimuth and/or 
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FIG. 15: Angular azimuthal correlations dcr/dcj)- at \fs = 1960 GeV for different (scalar) cuts As = 0, 1, 2, 3 GeV for NLO 
collinear (dashed), Kwieciriski (solid), BFKL (dashed), KL (dotted) and KMR (dash-dotted). Here pi,t,P2,t G (5,20) GeV and 
2/1,2/2 G (-5,5). 
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FIG. 16: Angular azimuthal correlations da/dcfi- at y/s = 1960 GeV for different (vector) cuts Ay = 0, 1, 2, 3 GeV for NLO 
collinear (dashed), Kwiecihski (solid), BFKL (dashed), KL (dotted) and KMR (dash-dotted). Here pi,t,P2,t G (5,20) GeV and 
2/1,2/2 G (-5,5). 



13 



in the two-dimensional (pi.t,P2,t) space which cannot 
be populated in the standard next-to-leading order ap- 
proach. In the latter case the fc t -factorization seems to 
be a useful and efficient tool. 

We believe that the photon-jet correlations can be mea- 
sured at Tevatron. At RHIC one can measure jet-hadron 
correlations for rather not too high transverse momenta 
of the trigger photon and of the associated hadron. This 
is precisely the semihard region discussed here. In this 
case the theoretical calculations would require inclusion 
of the fragmentation process. This can be done easily as- 
suming independent parton fragmentation method using 
fragmentation functions extracted from e + e~ collisions. 
This will be a subject of the following analysis. 
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B. Matrix elements for 2^3 processes 

In order to obtain parton-parton — > 7-jet matrix ele- 
ments for the next-to-leading order one can use the fol- 
lowing expression 

spins ^ col 

o(r< 1 at \ _■_ at °2«7 + a^a& 

2 (^F - TrA c )a 4 + N c 

2 a 9 

a\ + fly 



a 2 a 3 aQa 7 01030507 01020506 



for the 7(pi) + l{Pi) — * g{fa) + g(fa) + q(fa) process 
obtained in Q. Here Cf = 4/3, Nc = 3, where 



9l 

gL 



4na s (pl t ) 



and 



V. APPENDIX 

A. Matrix elements for 2^2 processes with initial 
off-shell partons 

In this paper we include four 2^2 processes such 
as qq -> jg,qq -> 75,50 -> jq,qg -> 70 important 
at midrapidity and relatively small transverse momenta. 
The corresponding matrix elements for the on-shell initial 
partons read 



7rn! emA /ai, a Q!2 t g(167r) (-1 ( | + , 



I X 99-79 1 2 



I -^9-9-79 1 2 



"<V-emV a l,a° t %s(l&K) ( S ) ( ~ + J ) , 



7ra em ^/ai iS a2, s (167r) 



9/ \u t, 
1\ (u s 

3 U + w ' 



\M q 



,„ \ <> 1 ,n 2 s (167r) ( -- J (j + 1 



The matrix elements for the off-shell initial partons were 
derived in Ref.Q- To a good approximation the ma- 
trix elements for the off-shell initial partons can be also 
obtained by using the on-shell formulae (fT8|) but with 
s,t,u calculated including off-shell initial kinematics. In 
this case s + i + u = k\ + k\, where k\,k\ < denote 
virtualities of initial partons. Our prescription can be 
treated as a smooth analytic continuation of the on-shell 
formula off mass shell. With our choice of initial parton 
four-momenta k\ — —k\ t and k\ — —k\ t . 

Explicit formulae for exact off-shell matrix elements 
were calculated and can be found in Ref. 0] • In this paper 
we compare results obtained with both (approximate and 
exact) ways. 



«i =P2 -Pi, a-5 

a 2 = P2 ■ fa, a e 

a* = P2 ■ fa, a 7 
a A =p 2 - fa, 



fa- Pi, a 8 = fa-pi, a 1Q = k 1 -p 1 , 
k 2 ■ ki , a 9 =fa - fa, 
k 2 ■ k 3 , 



are redundant invariants. The longitudinal momentum 
fractions are calculated as: 

xi = (pi,*e yi + P2 , t e V2 +p 3 ,te V3 )/V^ 
xi = (pi,te- yi + P2 , t e- V2 +P3, t e- V3 )/vG 

As an example the expression for matrix elements for 
the second diagram 

g(Pi) + g(P2) -> lifa) + q(fa) + q(k 2 ) 



in Fig 2(b) we get from 



liPl) + q(j>2) -> g{fa) + g(k^+q(k 2 ) 

diagram (see FigJTT]) if we make the following replace- 
ment 

Pi^ fa, 

P2 fa, 
fa -> P2, 

thus obtaining: 

a\ — > 09, 05 — > 06, a§ — > 02, aio — + aio, 

02 — ► as, ag — > 05, ag — ► ax, 

a 3 —> a 3 , a 7 —> a 4 , 
a 4 — > a 7 , 
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FIG. 17: Diagram of the yq — » ggg process. 
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